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Abstract 

For reaction-diffusion processes without exclusion, in which the particles 
can exist in the same site of a one-dimensional lattice, we study all the 
integrable models which can be obtained by imposing a boundary con- 
dition on the master equation of the asymmetric diffusion process. The 
annihilation process is also added. The Bethe ansatz solution and the 
exact iV-particle conditional probabilities are obtained. 



1 Introduction 

There are a variety of phenomena that can be explained by stochastic mod- 
els, and their non-equilibrium behaviors can be understood by rather simple 
rules [1-3]. One of the important examples of these models are the reaction- 
diffusion processes on a one-dimensional lattice for which their dynamics are 
fully specified by their master equation [4,5]. A simple example of reaction- 
diffusion process is asymmetric simple exclusion process (ASEP) [2,6,7], which 
is known to be relevant to various fields of science like the kinetics of biopolymer- 
ization [8,9], traffic models [10], polymers in random media, dynamical models 
of interface growth [11,12], noisy Burgers equation [13], study of shocks [14,15], 
sequence alignment [16], and molecular motors [17]. For recent reviews, see for 
instance [18-20]. 

The ASEP is a lattice model in which each particle hops to its right (left) 
nearest-neighboring site with a probability Dndt(DLdt) in an infinitesimal in- 
terval dt. In addition, particles are subject to hard-core exclusion: each site is 
either occupied only by one particle or empty. ASEP has been studied in [21] 
by introducing a master equation which describes the time evolution of proba- 
bilities P(x\, ■ ■ ■ t xpf;t), when the particles are not in neighboring sites, and a 
so-called boundary condition, which specifies the situations in which the prob- 
abilities go outside the physical region x\ < X2 < ••• < xn- These happen 
when some of the particles are in adjacent sites and the master equation can 
not be applied to them. It has been shown that the model is integrable in the 
sense that the iV-particle S'-matrix is factorized into a product of two-particle 
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S'-matrices. The coordinate Bethe ansatz has been used in this proof. Note that 
the S-matrix can completely determine the dynamics of a Markovian process, 
i.e. the iV-particle probabilities of a model. 

By choosing other suitable boundary conditions, without changing the mas- 
ter equation, one may study the more complicated reaction-diffusion processes, 
even with long-range interaction. In ref. [22], the so-called drop-push model has 
been studied by this method. In this model the particle hops to the next right 
site even it is occupied. The particle hops to this site by pushing all the neigh- 
boring particles to their next right sites, with a rate depending on the number 
of right neighboring particles. The generalization of this model, by consider- 
ing both the right and left hopping, has been done in ref. [23]. This method 
has been also applied to more-than-one-species situations, which become more 
complicated. The complexity arises from the above mentioned factorization 
of A-particle scattering matrix. In these cases, the factorization demands the 
two-particle S'-matrices to satisfy the spectral Yang-Baxter equation. Various 
solvable multi-species models have been studied in this way, of which the most 
recent general cases have been discussed in [24] and [25]. 

All of the above studies have been restricted to interactions which include 
the hard-core exclusion. This made some simplification. In ref. [26,27], an 
asymmetric diffusion model without exclusion has been shown to be integrable 
and to have the same R- matrix as that of the ASEP. In [28] , it has been shown 
that the processes: 

mn — > m — l,n + I with rate -rrr : 

r\ 

D R <y l 

mn — ► m-\-l, n— I with rate — t-. — , (1) 

r\ 

is integrable, in the sense of the above mentioned two-particle factorization. The 
numbers "to" , "n" , • • • indicate the particle numbers on a site, 7 = Dl/Dr, and 

1 - V 

[l] = ~j • (2) 

1-7 

Note that in [28], the reaction rates of eq.lQJ are scaled by D R . These processes 
are obtained by imposing the boundary condition 

P N (--- ,Xj,Xj-l,--- ;t) = D R P N (--- ,Xj-l,Xj-l,-" ;t)+D L P N (- ■ ■ ,Xj,Xj,--- ;t) (j = 1, ■•■ ,N-1), 

(3) 

on master equation of ASEP: 
d N 

—P(xi, ■ ■■ ,Xj,-" ,x N ;t) = ^D B P(x 1 , ■ ■ ■ ,Xj-!,Xj - l,Xj + i ■ ■ ■ ,x N ;t) 

N 

+ J~]DlP(xi,--' ,Xj-i,Xj + 1, atj+i • • • ,x N ;t) 
3=1 

-NP(xi, ■■■ ,Xj,--- ,XN]t), (4) 

in which we have used a time scale so that Dr + Di, = 1. In above equation, 
P(xi,--- , Xjv;i) is the probability for finding at time t the particles at sites 
X\, ■ ■ ■ ,xn- We take these functions to define probabilities only in the physical 
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region x\ < x% < • • • < xjv. In fact, in the domain f2jv = a?i < X2 < • ■ • < 
ijv C Z^, the function P is the probability defined above, whereas in Z N 
it is defined by the master equation Q , but it is not a probability. The master 
equation J3J is only valid for 

< Zj+i, (5) 

since for xi = aji+i, there will be terms with a;j = Xj+i + 1 in the right-hand 
side of eq.||3J|, which is out of the physical region. One can, however, assume 
that ||2J is valid for all physical region xi < Xi+i, and impose certain boundary 
condition for Xi — Xi+i. The boundary condition J2J leads to interactions lfT)l. 

In this paper we want to study all possible boundary conditions for single- 
species systems and derive all the integrable one-dimensional reaction-diffusion 
processes without exclusion which can be obtained by this method. The scheme 
of the paper is the following. In section 2, we show that there are two types 
of boundary conditions when Dr ^ and Dl 0. The first is one considered 
in [28], i.e. eq.©, which we call it the type 1 model, and the second one is: 

D R P(x, x-1) + D l P{x + 1, x) = P(x, x). (6) 

Here we have suppressed all the other coordinates for simplicity. We show that 
the interactions of this type 2 model are 

mn —> m — 1, n + 1 with rate Dr, 

mn —> m+l,n— 1 with rate D^. (7) 

This is the reactions restricted to I = 1. In other words, to have an integrable 
model, it is not necessary to have the simultaneous hoppings of any number of 
particles from a common site to the neighboring site ( as indicated in Q), but 
the one particle hoppings can also lead to integrable models. Note that the 
reactions (JJJ) is not a subset of processes 

In the totally asymmetric case with Dl — 0, we show in section 3 that 
there is a new boundary condition which is a linear combination of boundary 
conditions J2J and ©, i.e. 

P(x, x-l) = lP(x, x) + fiP(x -l,x- 1). (8) 

It is shown that this type 3 model contains the reactions 

mn —> m — l,n + 1, (9) 

with rates 

n = 1 : — ■ (io) 

! + £ + ■••+(' 

This is an interesting one-parameter family of interactions. 

In section 4, we generalize the boundary conditions ©, © and © to include 
annihilation to the processes Q), and respectively. The resulting models 
are rather involved. To be specific, we consider Dl — case of Q, and show 
that the boundary condition 

P(x,x- 1) = fiP(x- l,x- 1) (/i<l) (11) 
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describes the reactions 



mn — > m — /, n + I with rate /i , 
f n-l 

n — > < n — 2 with total rates n — Y^/A (12) 

• '=° 

The second reactions are the annihilation processes. We call this model as the 
type 4 model. It must be mentioned that we can not extend our investigation 
to include the creation processes. The main reason is that if we do so, the 
evolution equation of n-particle sector will not become closed and will depend 
on the more-than-n-particle configurations. 

The Bethe ansatz solution for different models is discussed in section 5 and 
the exact iV-point conditional probabilities of type 4 model is obtained in section 
6. Some interesting physical quantities are also obtained. Finally we discuss 
the multi-species extension of these processes in the last section and show that 
this generalization is not possible for the reactions without exclusion. 



2 Reactions with Di ^ and Dr ^ 

Consider the master equation Q for two particle sector, when the particles are 
at site x\ = X2 = x: 

P(x, x) = D R P(x-l, x)+D R P(x, x-l)+D L P(x+l, x)+D L P{x, x+l)-2P{x, x). 

(13) 

The second and third probabilities in the right-hand side of the above equation 
are out of the physical region, and must be defined through some boundary 
conditions. There are only two possibilities which are consistent in more-than- 
two-particle sectors. The first one is 

fi'P(x, x - 1) = D R P(x - 1, x - 1) + D L P{x, x), (14) 

and the second one is 

D R P{x, x - 1) + D L P{x + 1, x) = l'P(x, x). (15) 

In the first choice, we take any unphysical terms of eci. (|13J) as a linear combina- 
tion of physical functions, and in the second choice, we take the whole unphysical 
terms as a linear combination of the physical probabilities. The right-hand sides 
of eas. (|14[) and IjlfrJ) are the only allowed combinations which one can write. This 
was discussed in [25] for ASEP cases. In fact, one can use any other parameters 
a and (3 instead of D R and Dl in ea. 1)14(1 . with condition a + (3 = 1. Also one 
can add two other terms fiP(x — 1, x — 1) and vP{x + 1, x+ 1) to the right-hand 
side of eci. (|15f) . but it can be shown that for obtaining a consistent description in 
more-than-two-particle sectors, one must take [i = v = ( see section 3 of [25] 
for more details in ASEP case). 

To obtain the range of parameters // and I', one can use eq.(@J to show that 

d0dt^2 ^2 P(x 1 ,x 2 ;t) = -^P{x,x;t)+D R ^P(x,x-l;t)+D L ^P(x + l,x;t) 

X 2 X\<X 2 XX X 
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= -5^P(x,a?;t)+5^P(as ) a;-l;t). (16) 

X X 

For boundary condition l|14|l . ca. (|16|l results 

d0dtY^ P(.xi,x 2 ;t) = {l0fi' -l)^2P(x,x;t), (17) 

X2 Xi<l2 x 

and for boundary condition 1|15|) . it results 

c*06»t^ ^ P^!,^;*) = (Z'-l)^P(z,z;t). (18) 

X2 Xi<X2 X 

In the first step, let us exclude the annihilation processes and therefore it is 
assumed that the number of particles is constant in time. In this way, the 
eqs.JITJ) and ^HJ lead to 

fjf = 1 (for boundary condition (|14l> ). 
I' = 1 (for boundary condition (|15ll ). (19) 

Ea. (|14fl exactly becomes the one studied in [28], i.e. eq.Q, and it induces 
the reactions but the second one, ea. (|15fl . is new. To obtain the reaction 
introduced by this boundary condition, with V — 1, we first consider P(x,x) in 
ea. (fH5|l . Using cas. ((T5|) and lfT^|l . one finds 

P(x, x) = D R P{x - 1, x) + D L P{x, x + 1) - P(x, x), (20) 

which is the evolution equation of the following two-particle reactions 

10 -> 01 with rate D R , 

01 10 with rate D L . (21) 

Remember Dr + Dl = 1- Generally, for n particles existing at a common site 
x, one finds from |0J: 

d 

— P( x, ■ ■■ ,x ) = DrP(x - 1, x, ■ , x ) + DlP(x, ■■■ ,x,x + l) 

n 71 — 1 71—1 

n 

+ DrP(xi = X, • • • , Xj— l — X^Xj = X — 1, = x • • • , x n = .x) 

J=2 

+ DlP(xi = X, ■ ■ ■ , = X, Xj-i = X + 1, Xj = X, ■ ■ ■ , x n = x) 

i=2 

—nP(x, ■ ■ ■ ,x) 

n 

= D R P(x-l,x,--- ,x) + D L P(x,--- ,x,x + l)-P(x,--- ,x), (22) 

n — 1 n — 1 

which is obviously the evolution equation of reactions (JjJ) . 

3 Totally asymmetric diffusion with Dl = 

In Dl = 0, the boundary condition © and © becomes 
P(x,z:-1) = P{x - l,x - 1), 
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P(x,x-1) = P{x,x), (23) 
respectively, which describe the interactions 

mn — > to — l,n + I with rate 1 (24) 

and 

mn — > to — 1, n + 1 with rate 1, (25) 

respectively. But as was first noted in [22] for ASEP, the linear combination 
of these two boundary conditions may result a new integrable model. So we 
consider 

P(x, x-l) = lP(x, x) + [iP(x - l,x - 1), (26) 
as the boundary condition for the following master equation 

d N 

—P(xi, ■ ■ ■ ,x N ;t) = ^2 [P(xi, ■ ■ ■ ,x 3 -i,Xj - l,x j+ i ■ ■ ■ ,x N :t) - P(xi, ■ ■ ■ ,x N ;t)] . 

(27) 

First we note that 

d0dtJ2 P{xi,x 2 ;t) = {l + ii-l)^2P{x,x;t), (28) 

X2 X\<X2 X 

from which the conservation of number of particles results 

l + H = l. (29) 

To obtain the resulting reactions, besides the diffusion 10 — > 01 with rate 1, we 
first consider P(x, x). Using l|2t)|l and (|27|l . it is found 

P(x, x) = P(x -l,x) + {iP{x - 1, x - 1) - (2 - l)P(x, x), (30) 

which represents the reaction 20 — > 02, with rate fi = 1 — I, as the source and 
sink of this state. To find the reactions in general case, we first prove a lemma. 
Lemma: Equation l|26|l implies, for arbitrary n, the following 

P( x, ■ • • , X , x - 1) = r n+1 P{ x- ,x- 1 ) + (1 - r n+ i)P( x, ■ ■■ (31) 

n n+1 n+1 

with r n defined through l|10|l . 

Proof: We proceed by induction. For n = 1, (|31|l reduces to l|2t)|) . Assuming 
(131 j l is correct for n — 1, then using l|26|l . we have 

P(x, ■ ■ ■ , x, x — 1) = lP(x, ■ ■ ■ ,x) + ^P{x, ■ ■ ■ , X, x — 1, x — 1) 

n+1 — 1 

= ZP Qc, ■ ■ ■ ,x ) + fj.[r n P( x ■ ,x - 1 ) + (1 - r n )P( x, ■ • • - 1)], 

n+l n+1 n 

(32) 



so 



P(x, ■ ■ • ,x,x - 1) = r n+1 P( x ,x-l ) + s n+1 P(x, ■ ■■ ,x), (33) 

n+l n+l 
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where 

r n+ i = /xr„0l - - r n ) , s n+ i = l0l - - r„). (34) 

It is seen that r n +i + s n +i — !• The value of r n can be found by solving the 
first equation of 134(1 . which can be written as 

r~li = l + Wr- 1 . (35) 

Using r^ 1 = l0fx = 1 + l0fi, (|3*5)) leads to ea . (fTT^I for r„s. This proves the 
lemma. ■ 

We now consider the evolution of P(x, ■ ■ ■ ,x). Using (|2*7|) and (31), we find 



n 

P(x, ■ ■ ■ , x) = P(x\ = x, ■ ■ ■ , Xj-i = x, Xj — x — 1, Xj+i — x ■ ■ ■ , x n — x) — nP(x, ■ ■ ■ ,x) 



Of 



= ^^[rjP(xi = x — 1, • ■ • , Xj = x — 1, xj + i = x ■ ■ ■ ,x n = x) + (1 — rj)P(x, • • • , x)} 

n 

—nP(x, ■ ■ ■ ,x) 

n 

n n 

= rjP(xi = x — 1, ■ • • ,Xj = x — 1, a^j+i = or • • • , a;„ = x) — (V^ rj)P(x, ■ ■ ■ ,x). 
i=i j'=i V v 

(36) 

This equation shows that a collection of j particles hops from site x — 1 to 
site x with rate r*j , and this proves that the induced reactions of the boundary 
condition (8) are those indicated in (9). 

At I = 1 (fi = 0), eq.(10) results r; = 5n. So the reactions are 

ran — » m — 1, n + 1, with rate 1. (37) 

This is nothing but eq.(25). At I = Q (/i = 1), eq.(10) gives r; = 1 for all 
Zs. So the reactions are those shown in (24). In fact, at I = 0, all the multi- 
particle hoppings occur with equal rate 1. By increasing I, the greater number 
of simultaneous hoppings happen with lower rates, until at I = 1, in which only 
the one-particle hopping is allowed. 



4 Annihilation-diffusion processes 

Adding annihilation to the previous reactions, results in the decreasing of the 
number of particles with time. Therefore the parameters //, V and (/, /i) in 
eqs.(14), (15) and (26) become 

// > 1 (for boundary condition (|14|l '). 
I' < 1 (for boundary condition Q\lb]i ). 
fi + 1 < 1 (for boundary condition (|2*6^l ). (38) 

These are obtained from eqs.(i7), (18) and (28), respectively. Note that the 
annihilation interactions only appear in sink terms of the master equation, since 
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if we consider an initial state with n particles, no annihilation process can lead 
to a n-particle state at any later time t. 

To obtain the corresponding interactions of either of boundary conditions, 
we begin with (15) with I' < 1. Considering P(x,x), one finds 

P(x, x) = D R P(x -1,x) + D l P(x, x + 1) - (2 - l')P(x, x), (39) 

which shows the annihilation rate 1 — V for the state of two particle at a common 
site. In general case, one finds, similar to eq.(22): 

d 

— Pix^^x) = D R P(x-l,x 1 -^x)+D L P(xr^ ^x,x+l)-[n-(n-l)l']P(x, - ■ ■ ,x), 

n n—1 n—1 

which is the evolution equation of the following reactions: 

ran — ► m— l,n+l with rate Dr 1 
mn — ► m + l,n— 1 with rate Dl, 
n—1 



(40) 



n 



n - 2 with total rates (n-l)(l-Z')- (41) 



In the case of boundary condition (26) with I + [i < 1, it is not easy to obtain 
a compact form for interaction rates in general case. For example in two- and 
three-particle sectors, considering P(x,x) and P(x,x,x), one finds 

mn — > m — l,n+l with rate 1, 

ran — > m — 2, n + 2 with rate /z, 

mn — > m — 3,n + 3 with rate a 2 0l — la 1 

2 — > 1 with rate a = 1 — (I + fi), 

2 

~ with total rates a(u 2 + ol\i + 2)0l — Z/i. (42) 

In this case, we have a two-parameter family of interactions. 

When we consider the boundary condition (14) with // > 1, again we can 
not find the compact relations and it is better to restrict ourselves to the subset 
of totally asymmetric reactions with Dl — 0. This means that we take 

P(a;,a;-l) = /iP(ar-l,a;-l) (/x < 1) (43) 

as the boundary condition ( u = 1/// ) . Considering the general case P(x, ■ ■ ■ ,x), 
results 



d 

( X . ' ' ' ■ X^j — ^ ^P{^X~y — X j ' ' ' j Xj ~ 1 — X ^ X j — X ' 1^ — X ' ' ' "i Xn — X^j Tl^P^X^ ' ' ' . X^j 

V 7 — 1 

n J n 

n 

= fjp~ 1 P(x\ = x — 1, • • • , Xj = x — 1, = x • • • , x n = x) — nP(x, • • • jx). 

(44) 

This equation shows that the sources are the simultaneous hoppings of j particles 
(j = 1, • • • , n) from the common site x — 1 to x, with rates \i 3 ~ . So the sinks 
are also these hoppings and the remaining rate of the sink terms in eq.(44), i.e. 
n — ^ i IS due t° t ne annihilations of these n particles at site x. So, it is 

proved that the boundary condition (43) induces the reactions (12). 
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5 Bethe ansatz solution 

Now we try to solve the resulting evolution equations, in all the discussed cases, 
by the following Bethe ansatz 

P(x;i) = e- B "VW, (45) 

where x = (x\, • • • , Xn), and 

1>(x) = J2 A * eia{p) '*- ( 46 ) 

The summation runs over the elements of permutation group of N object. In- 
serting (46) in master equation (4), results 

N 

e n = E C 1 - D * e ~ vPk - DLeWk ) ■ ( 4? ) 

fe=l 

To determine A a , we must insert the expression (46) in the boundary conditions. 
For type 2 model, for instance, we must use eq.(6), which results 

D R tp(- ■ ■ ,Xi = x, x i+ i =x—l,--- )+D L ip(- ■■ ,Xi = x+1, Xi+i = x,- ■ ■) = ip(- ■ ■ ,Xi = x, x i+ i = x, ■ ■ ■ ), 

(48) 

and using (46), gives 



A a a k — 0. 



DRe -ia( Pk+1 ) + Dhe i<vu) -i A a + D R e- ia{Pk) + D ie lCT (P*+i) - 1 

(49) 

(Tfe is an element of permutation group which only interchanges pt and Pk+i, 

°k ■ (pi, • • • ,Pk,Pk+i, ■■■ ,Pn) -» (pi, ■ ■ ■ ,Pk+i,Pk, ■ ■ ■ ,Pn)- (50) 
Eq.(49) gives A aak in terms of A a as following 

A aak =S^(a(p k ),a(p k+1 ))A a (51) 

where 

S^\z u z 2 ) = -D R z 2 + D L z^ - l0D RZl +D L z^ - 1, (52) 

where z k = e~ lPk . Eq.(51) allows one to compute all A a s in terms of A\, which 
is set to unity. 

The same procedure can be applied to other boundary conditions. For ex- 
ample for type 3 model, boundary condition (8), one finds 

S {3,) (z x ,z 2 ) = -l + n(zx + z 2 ) - z 2 0l + n(zi + z 2 ) - zi, (53) 

and for boundary condition (15), with I' < 1, 

S( Zl ,z 2 ) = -D R z 2 + Dlz^ 1 - l'0D RZl + Dlz^ 1 - I'. (54) 

These solutions can be used, in principle, to calculate the conditional probabili- 
ties P(xi, ■ ■ ■ ,xn] t\y\, • • • , un; 0). This is the probability of finding the particles 
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at time t at sites x\, ■ ■ ■ , xn, if at t = 0, they were at sites j/i, • ■ ■ , Un, respec- 
tively. But unfortunately, the standard method, used for example in [21-25], 
can not be used here. This is because the initial condition 

P(x;0|y;0) = <5 x , y , (55) 

is satisfied by standard expression 

P(x; t\y; 0) = J d N P 0{2n) N e - Bw 'e- ip - y V(x), (56) 

only when t/j < yi+i and Xi < Xi+i. This condition only satisfied by exclusion 
processes and for processes without exclusion, we must look for other methods. 



6 iV-particle conditional probabilities 

In some special cases, it is possible to calculate the conditional probabilities in 
terms of a specific determinant. This was first proved in [21], and then used 
for other cases in [22,23]. Now we want to check that this method does work 
here, and as a specific example, we consider the type 4 model with the boundary 
condition (11). 

We set the following ansatz for TV-particle conditional probabilities 

P(x; t\y; 0) = e- Nt det[G(x; t\y; 0)], (57) 

where G is a N x N matrix with elements 

G y (x;t|y;0) = gi-j(xi - yf,t). (58) 

Inserting (57) in eqs.(27) and (11), results 

d0dtG t (x;t) = Gi(x-ht), 

Gi- X {x;t) - iiGi-x{x-\;t)+pGi{x-\;t), (59) 

where Gi denotes the z-th row of matrix G, and (3 is an arbitrary constant. In 
terms of functions g p (n; t), (59) becomes 

d0dtg p (n;t) = g p (n — l;t), 

g p (n;t) = \ig p {n - 1; t) + /3g p+ i(n - 1; t). (60) 

Introducing the z-transform 

oo 

g p (z,t)= zn 9 P (n;t), (61) 

n— — oo 

we find 

d0dtg p (z,t) = zg p (z,t) => g p (z, t) = e zt g p (z, 0), 

g p+1 (z,t) = 1- fj,z0(3zg p (z,t). (62) 

The second equation yields 

g p {z,0) = (1- ^z0f3z) p g o (z,O). (63) 
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Using P(x\ 0\y; 0) = go(x — y; 0) = S X:V , one finds go(z, 0) = 1. Finally 

~g p (z 7 t) = e zt {l-nz0z) p , (64) 

in which we choose (3 — 1. The functions g p (n; t) can be obtained by expanding 
the generating functions g p (z, t). For p > 0, the expansion yields 

g p (n;t) = £ ( £ ) (^) m ^->+ P -m)!, (65) 

m=0 ^ ' 

and for negative p, it yields 

fl_ H (n;t)= £ ( W+ ™ )M ro *"- W - ro 0(n-N-m)!. (66) 

m=0 ^ ' 

We have thus obtained the explicit relation for conditional probabilities. 

It can be checked that the resulting function satisfies the desired initial 
condition. At t — 0, where Xi = yi, we have g_\ p \(Q;t) ~ and go(n;0) = <5„.o, 
which result P(x; 0|x; 0) = 1. Also it may be interesting to obtain the rate of 
decay of a delta function distribution. Suppose at t — 0, there are N particles 
at the same site y. We want to obtain the probability of finding all the particles 
at their initial positions at later time t, i.e. P(y; t\y; 0), where y = (y, ■ ■ ■ , y). 
Using (7_| p |(0;i) = and go(Q',t) = 1, we obtain 

P(y;t|y;0) = e~ m det[G(y;t|y;0)] 

= e- Nt , (67) 

which is independent of /z! It is an interesting result. At \x — 1, where there is no 
annihilation, the rate of simultaneous hoppings of particles are all 1, irrespective 
of the number of particles. By decreasing /x, the rates of simultaneous hoppings 
decrease with increasing the number of particles (see eq.(12)). At the same 
time, the rate of annihilation increases in such a way that the total decaying 
rate remains constant. 

The above determinant method can be applied to any model which its bound- 
ary condition equation contains only two terms. Otherwise it does not lead to 
the consistent relations for matrix elements g p (n\ £)'s. So for Dr = or Dl = 
cases of boundary conditions J2J and © , and A = or pi = cases of boundary 
condition JHJ, this method leads to explicit expressions for iV-particle condi- 
tional probabilities. 



7 Conclusion 

In the previous sections, we use all the allowed generalization of boundary condi- 
tion in the asymmetric reaction-diffusion processes without exclusion, to obtain 
several new integrable models for one-species cases. The further natural gener- 
alization , which has been discussed in many papers for the ASEP case, is the 
multi-species extensions of these models. In multi-species studies, one considers 
a p-species system with particles A\, - ■ ■ , A p . The basic objects are the proba- 
bilities P ai ... aN (x\, ■ ■ ■ ,XN;t) for finding at time t the particles of type ot\ at 
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x\, particle of type a 2 at x 2 , etc. The master equation, instead of (4), is [25] 

JV 

^ ] DnP ai ... aN (x\, ■ ■ ■ , Xj-i, Xj — 1, xj+i ■ ■ ■ , xn; t) 

N 

,Xj—i,Xj + 1, xj+i • • • , xn) t) 

-JVP ai ... ajv (a;i, • • • • • • ,x N ;t). (68) 

Now if we want to use this equation for reactions without exclusion, the problem 
will arise when some of the particles are in the same position. Consider, for 
example, the two-species case P aia2 (x, x). It is seen that eq.(68) has the term 
P ai a 2 { x ~ IjE) as the source term of the desired state, but does not contain 
Pa 2 a 1 {x — l,x), which is as important as the first term. In fact, this is the 
source of many difficulties arise in multi-species extension of reaction-diffusion 
processes which have not exclusion. So it seems that the integrable models 
discussed in the previous sections are all that one can obtain in this context. 

Acknowledgement: I would like to thank the research council of the Uni- 
versity of Tehran for partial financial support. 
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